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Abstract 

Deep inelastic charged-current reactions have been studied in e + p and e~p col- 
lisions at a center of mass energy of about 300 GeV in the kinematic region 
Q 2 >200GeV 2 and x>0.006 using the ZEUS detector at HERA. The integrated 
cross sections for Q 2 >200GeV 2 are found to be (T e + p -&x = 30.3 -2.6 P' 3 an< ^ 
(T e - P ^ u x = 54.7 1 |'| pb. Differential cross sections have been measured as 
functions of the variables x, y and Q 2 . From the measured differential cross 
sections da/dQ 2 , the W boson mass is determined to be My/ = 79^14 GeV. 
Measured jet rates and transverse energy profiles agree with model predictions. A 
search for charged-current interactions with a large rapidity gap yielded one can- 
didate event, corresponding to a cross section of a e +p^x{Q 2 >^^ GeV 2 ; rj max < 
2.5) = 0.8 to* ±0.1pb. 
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1 Introduction 



In comparison to fixed target neutrino scattering experiments |l[], the HERA ep collider 
extends the kinematic region for studying charged-current (CC) deep inelastic scattering 
(DIS) by about two orders of magnitude in Q 2 . In addition, HERA allows measurements 
of CC DIS at lower x (x < 0.01). Here Q 2 is the negative square of the four-momentum 
transferred between the electronS and the proton, and x is the Bjorken scaling variable. 

Both the HI and ZEUS || collaborations have previously reported cross section mea- 
surements for the process e~p — > u e X, where X denotes the hadronic final state. These 
investigations, based on the data collected in 1993, established that the Q 2 dependence 
of the CC cross section is consistent with that of the W propagator and that the CC 
and neutral-current (NC) cross sections have similar magnitude for Q 2 ^ [Mw de- 
notes the W mass). The HI collaboration has also measured the integrated cross section 
|4j and the differential cross section da/dQ 2 |J for e + p — ► V e X with missing transverse 
momentum (f t ) above 25 GeV. 

In 1994 ZEUS has collected 2.93 pb -1 of e + p data and 0.27 pb -1 of e~p data in collisions 
of 820 GeV protons with 27.5 GeV electrons. In order to study both e + and e~ induced 
CC reactions, these data samples have been combined with the 0.55 pb _1 of e~p data 
taken in 1993 with an electron beam energy of 26.7 GeV ||. 

In the quark parton model (including the leading order QCD evolution in Q 2 ), the unpo- 
larized e + p and e~p charged-current differential cross sections are:i 

2 _ 

£ [xu t (x,Q 2 ) + (l-y) 2 xd r (x,Q 2 )] , (2) 



rl 2 rr H 2 / M 2 

a u e~p^uX " p I 1V1 W \ \ ^ f / /o2n /i \2 



dxdQ 2 2vr \M^ + Q 2 

where Ui and di are the densities of up-type and down-type quarks of the i th generation 
in the proton, Gf is the Fermi coupling constant, and y is the fractional energy transfer to 
the proton in its rest system. The factor (1 — y) 2 , which suppresses the quark (antiquark) 
contribution to the e + p (e~p) cross section, is a direct consequence of the V— A structure 
of the weak coupling. 

As equations [I] and |2] show, the charged current couples to different quark flavors for e + 
and e~ beams and the valence and sea contributions have very different y dependences. In 
addition, the Q 2 dependence of the cross section allows a measurement of the W propaga- 
tor effect and thus a determination of the W mass in the space-like region, complementing 
the measurements of direct W production at pp colliders 0. 

This paper reports on measurements of the integrated and differential cross sections for 
both e + p and e~p CC DIS in the kinematic region Q 2 > 200 GeV 2 . The differential cross 
sections were measured as functions of x, y and Q 2 , where the latter dependence was 



1 In the following, "electron" is generically used to denote both electrons and positrons. 
2 For the current integrated luminosity, any contributions to the CC cross section from hard subpro- 
cesses involving 3 rd generation quarks (b, t) can be neglected. In this approximation the relation between 

quark densitites and CC structure functions is given by F% p ~^ vX oc d+s+u+c, F% v ^ vX oc u + c+d+s, 

xF% p^ vX (xd+s — u — c and xF^ p~* 1 ' x oc u + c — d — s. 
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used to determine the W mass. The available statistics allowed studies of the hadronic 
final state at high Q 2 in CC DIS: the jet rates and transverse energy jet profiles were 
determined and compared to Monte Carlo predictions. The number of events with two 
jets (plus the proton remnant) directly measures the rate of hard QCD subprocesses. 
Recently, diffractive NC DIS has received particular attention These reactions 

are characterized by an absence of final state particles in a wide rapidity interval between 
the outgoing proton system (which escapes through the forward beam hole) and the rest 
of the hadronic final state. Such a rapidity gap results in a quiet region in the forward part 
of the main detector. Diffractive CC reactions are interesting because they predominantly 
occur at high Q 2 and also because they provide the possibility to study the flavor content 
of the states exchanged between proton and W boson. A search for events with a rapidity 
gap in our CC samples yielded one candidate, which will be described in detail. The rate 
of large rapidity gap CC events was compared to model predictions and to the NC case. 
After a brief summary of the experimental setup in section 2, the kinematics and the 
selection of CC events are described in section 3. The cross section calculations and the 
evaluation of the systematic uncertainties are outlined in section 4. Finally, section 5 
contains the results, which are summarized in section 6. 



2 Experimental Setup 

In 1994, a total of 153 colliding bunches were stored in the electron and the proton 
beams of HERA, together with an additional 17 proton and 15 electron non-colliding 
bunches which were used to study beam induced backgrounds. The r.m.s. bunch length 
was about 20 cm for protons and 1 cm for electrons, resulting in an event vertex position 
distribution^ with an r.m.s. width of 12 cm in Z. The typical instantaneous luminosity was 
1.5 • 10 30 cirr" 2 s -1 . The following description refers to the 1994 running period. Details 
of the 1993 experimental setup and HERA running conditions can be found in [TT| . 



2.1 The ZEUS Detector 

A description of the ZEUS detector is available in [|T2|, [13] . The primary components used 



in this analysis were the uranium-scintillator calorimeter [[14] and the central tracking 
detectors. 

The calorimeter is divided into three parts, forward (FCAL) covering the polar angle@ 
interval 2.6°<#<37°, the barrel (BCAL, 37°<#<129°) and rear (RCAL, 129°<#<176.1°). 
The calorimeters are subdivided into towers which subtend solid angles between 0.006 
and 0.04 steradians. Each tower is longitudinally segmented into electromagnetic (EMC) 
and two hadronic (HAC) sections (one in RCAL). The electromagnetic layer of each 
tower is further subdivided transversely into four cells (two in RCAL). The total number 
of cells is 2172, 2592 and 1154 in FCAL, BCAL and RCAL, respectively. Under test 



beam conditions, the calorimeter has an energy resolution of <je/ E=18%/ J E(GeV) for 



electrons and oeI E=35%/ \J -E'(GeV) for hadrons. The time resolution is below Ins for 
energy deposits greater than 4.5 GeV. 



3 The ZEUS coordinate system is right-handed with the Z axis pointing in the proton direction, 
hereafter referred to as forward, and the X axis horizontal, pointing towards the center of HERA. 
4 The polar angle 8 is defined here with respect to the nominal interaction point. 
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Energy which penetrates through the CAL can be measured in the Backing Calorimeter 
(BAC) which consists of proportional chambers interleaved with the iron plates which 
form the return yoke of the solenoid [|15| B 

The tracking system consists of a vertex detector |TI| and a central tracking chamber 



(CTD) |nj operating in a 1.43 T magnetic field parallel to the beam axis. The polar angle 
region in which the CTD allows accurate momentum measurement is 15°<#<164°. The 
transverse momentum resolution for full length tracks is cr{pt)/pt= [0.005]^ (GeV)] ©0.016. 
The luminosity was measured by the rate of high energy photons from the process ep — > 



ep7 detected in a lead-scintillator calorimeter [18] located at Z = —107 m. 



2.2 Trigger Configuration 



Data were collected with a three-level trigger system fL3[ . The first-level trigger condi- 
tions used in this analysis were based on electromagnetic energy, transverse energy, and 
missing transverse energy in the calorimeter. The thresholds were significantly below the 
corresponding cuts applied in the subsequent selection. The second level trigger rejected 
backgrounds (mostly p-gas interactions and cosmic rays) for which the calorimeter timing 
was inconsistent with an ep interaction. Events were accepted as CC candidates if the 
following conditions were all fulfilled: 

1. The missing transverse momentum f exceeded 9 GeV (see equation [|). 

2. f evaluated from all calorimeter cells except those adjacent to the forward beam 
pipe exceeded 8 GeV. 

3. Either a track was found in the CTD or at least 10 GeV was deposited in the FCAL. 

The third-level trigger applied stricter timing cuts and also pattern recognition algorithms 
to reject cosmic rays. 



2.3 Monte Carlo Simulation 

Simulated CC DIS events with electroweak radiative corrections were generated using 
lepto [|19j interfaced to HERACLES [fMj via django pT|] . The MRSA set of NLO parton 



density parameterizations was used, which is based on preliminary 1993 F2 measurements 
by ZEUS and HI as well as on recent data on W asymmetries from CDF and on Drell-Yan 
cross sections (cf. |2!| and the references therein). The latter measurements constrain the 
difference u — d, to which CC reactions are particularly sensitive. The hadronic final 
state was simulated using the color-dipole model as implemented in ARIADNE |23j for the 



QCD cascade and jetset |24[] for the hadronization. In addition, MC samples using the 



meps option of lepto instead of ARIADNE were generated. Each of the CC MC samples 
contained about 10 4 events. The results were corrected for radiative effects using these 
CC MC sets! 

MC samples of NC DIS and photoproduction events were employed to simulate back- 
grounds to CC DIS. For the NC DIS simulation the same programs were used as for 



5 Due to the present limited understanding of the Monte Carlo simulation of energy deposits in the 
BAC, this component has been used only for systematic checks in this analysis. 

6 All kinematic quantities are inferred from the four-momentum of the propagator. 
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the CC events. Samples of both direct and resolved photoproduction events were gener- 
ated using both pythia [[24] and herwig [25]. Photoproduction of cc and bb pairs was 



simulated using both pythia and aroma [E6[ . 

To calculate efficiencies for diffractive CC interactions the pompyt 
program was used. This program is based on the factorizing model of Ingelman and Schlein 



271 Monte Carlo 



28] and assumes a hard quark density in the pomeron which satisfies the momentum sum 



rule. 

All simulated events were passed through a geant based |29|] detector simulation and 
processed with the same analysis programs as the data. 



3 Kinematics and Event Selection 

In CC ep reactions the final state neutrino remains undetected. This fact imposed special 
conditions both for the reconstruction of kinematic variables and for the event selection. 

3.1 Kinematic Variables 

The kinematics of the interaction ep — > vX are defined by the four-momenta k and P 
of the incident electron and proton respectively, and the four-momentum of the final 
state neutrino (k') or the final state hadronic system (P'). The four-momentum transfer 
between the electron and the proton is given by q = k — k' = P' — P. 
In addition to s = (k + P) 2 , the square of the ep center of mass energy, three Lorentz 
invariants are defined in terms of these four-momenta: 

• Q 2 = —q 2 , the negative square of the four-momentum transfer, 

Q 2 

• x = , the Biorken scaling variable, 

2q ■ P ' J & 

q-P 

• y = , the fractional energy transfer to the proton in its rest system. 

k • P 

Neglecting mass terms, these variables are related by Q 2 = xys. 

The Jacquet-Blondel method |Kj was used to reconstruct the kinematic variables from 
the measured missing transverse momentum, 



\ 



(3) 



and the quantity 5, given by 

* = Pz) , (4) 



where the sums run over all EMC (HAC) cells with energy deposits E l above 60 MeV (110 
MeV) and the p 1 are the momenta assigned to each calorimeter cell (calculated assuming 
zero mass and using the geometrical cell center and the measured vertex position). The 
Jacquet-Blondel estimators of y, Q 2 and x are given by 

^ D2 ft 2 - P 2 fK\ 

2£ e I -Vjb syj B {l-yjB> 

where E e is the electron beam energy. 
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3.2 Selection of Charged— Current Events 

The number of CC triggers was about 10 5 , almost entirely due to p-gas interactions, 
cosmic rays, and beam halo muons. After a preselection which required a tracking vertex 
and stricter timing cuts than at the trigger level, about 6300 e + p and 1500 e~p candidates 
remained (cf. table p. In the following, we describe the cuts which were applied to extract 
the CC signal. Table [I] gives the numbers of events remaining after each selection cut, as 
well as the corresponding efficiencies estimated using the ARIADNE MC samples. 
The following conditions were imposed on all events passing the trigger: 

• f > 1 1 GeV was required. The trigger simulation indicated that this cut was 
sufficiently far above the 9 GeV trigger threshold to ensure high trigger efficiency. 

• The events had to have a tracking vertex with | Z vtx | < 45 cm, where Z ytx = 
at the nominal interaction point. This cut eliminated a large fraction of non-ep 
background. 

• ft° ui ' I 'ft > 0.6 was required, where f ont is the net transverse momentum for all 
calorimeter cells with a polar angle above 9° with respect to the nominal vertex 
position. This cut rejected p-gas and p-beampipe collisions, for which f is concen- 
trated at small polar angles. 

• There had to be at least one reconstructed track originating from the vertex, which 
had a polar angle between 15° and 164° and a transverse momentum exceeding 
0.2 GeV. This requirement removed cosmic ray events and p-gas interactions with 
spurious vertices caused by low-energy secondary interactions in the beam pipe. 

• The difference A<p between the azimuth of the net transverse momentum as mea- 
sured by the CTD tracks with polar angle between 15° and 164°, and the azimuth 
measured by the calorimeter, was required to fulfill | A0| < 1 rad. This requirement 
removed overlays of cosmic rays on ep-interactions. 

• In addition, a pattern recognition algorithm based on the topology of the calorimeter 
energy distribution was applied to reject cosmic rays and beam halo muons. 

Simulation of backgrounds due to NC DIS and photoproduction interactions showed that 
such events passing the above cuts were concentrated at low f t . The same was found for 
the non-ep background. Hence the following additional cuts were applied to events with 
f < 30 GeV: 

• Vjb < 0.8, which reduced NC DIS background. 

• ft/ E t > 0.4, where E t = Y,i \J (Px) 2 + (Py) 2 i s the total transverse energy. This cut 
demanded an azimuthally collimated energy flow and rejected photoproduction and 
also p-gas background. 

• If™** /ft > o.l, where P* racks = J (j2j P£ ackj ) 2 + (j2j P^ racki ) 2 and j runs over 
all vertex-fitted tracks with polar angle between 15° and 164°. This cut effectively 
tightened the ft° ut /ft requirement (cf. above) and also removed events with addi- 
tional non-ep related energy deposits in the calorimeter (mainly cosmic rays). 
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The resulting sample was visually scanned and 3 non-ep background events were identified 
in the e + p sample and removed (all were overlays of muons and p-gas reactions). No non- 
ep background was found in the e~p sample. The distribution in x and Q 2 of the final 
samples of 56 e + p events and 30 e~p events is shown in figure [l]. 



3.3 The Neutral— Current Control Sample 

The acceptance corrections and the methods of reconstructing the kinematic variables 
were based on the MC simulation described in section 2.3. In order to verify the accuracy 
of the simulation, several systematic checks were performed using a sample^ of high-Q 2 
NC events selected as follows: 

• As in the CC case (cf. section 2.2), the NC trigger requirements were based on 
the calorimeter energy deposits. The same timing cuts were applied as for the 
CC trigger. Photoproduction background was efficiently suppressed by the cut 
5 > 25GeV (cf. eq. |). 

• In the offline NC selection, an electron candidate with at least 10 GeV energy had 
to be identified by a neural-network algorithm using the pattern of the calorimeter 
energy deposits f31fl . The electron had to be isolated (less than 2 GeV of energy 
not associated with the electron in an (77, 0)-cone! of radius R = 0.5 centered on 
the electron direction) and had to have a matching track if it was in the region 
20° < 9 < 160°. In addition, the event had to have a reconstructed tracking vertex, 
and satisfy 35 GeV < 5 < 65 GeV and Q 2 > 200 GeV 2 , where Q 2 was reconstructed 
using the double angle (DA) method |32| . 



• For the selected NC events, the track(s) and the calorimeter deposits associated 
with the electron were deleted from the data, and the reconstruction of the event 
vertex was repeated. 

• These modified NC events were then processed through the complete CC analysis 
chain. In order to account for the different cross sections, a weight 

da e + p _ vx /dxdQ 2 , da e - p ^ uX /dx dQ 2 

w = 1 Ta a^2 for e P res P- w = ~a Ta ~Ara for e P 

da e + p ^ e +x/dx dQ z da e - p ^ e -x/dx dQ z 

was assigned to each event, where leading order differential cross sections, evaluated 
at x and Q 2 as reconstructed by the double angle method, were used. 

After the CC selection, 2269 e + p and 167 e~p events with Q 2 > 200 GeV 2 remained, 
where Q 2 was reconstructed as described in section 4.1. The respective sums of 
weights were 50.8 and 6.2. 

The control sample was not corrected for the efficiency of the NC event selection, for 
remaining backgrounds, or for migrations due to mismeasurement of x and Q 2 . The 
distributions of Q 2 B , xjb, Vjb and f t agreed sufficiently well with the corresponding 
CC MC distributions to allow for quantitative checks of the efficiency calculations and 
correction methods. 



7 Data collected in 1994 were used. 

8 The pseudorapidity r\ is defined as — ln[tan(6>/2)], where 9 is the polar angle with respect to the 
proton beam direction, taken from the reconstructed interaction point. 
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4 Cross Section Evaluation 



The integrated and differential cross sections were calculated from the observed numbers 
of events and the integrated luminosity after subtracting the estimated photoproduction 
background and correcting for detection efficiencies and migrations. 
The e~p data from 1993 and 1994 were combined by appropriately weighting the 1993 
event numbers in each individual bin under consideration. The weights, which account 
for slight differences in the electron beam energies and in the trigger and detector config- 
urations, were in the range between 0.95 and 1.05. 

All statistical errors have been evaluated using the observed numbers of events and asym- 
metric Poisson confidence intervals. 



4.1 Reconstruction of Kinematic Variables 

Due to the energy loss of the hadrons in inactive material in the central detector, the 
calorimeter measurement underestimated both f and 5. The detector simulation was used 
to derive corrected values ft cor and y COI , as second-order polynomials of the respective 
raw values f and yjs, where the coefficients for the y correction depend on f t . The 
correction of f t was about 20% at low f t , decreasing to ~ 10% for f t > 50 GeV; in y, the 
corrections were highest for y w 0.1 and low f (~ 20%), decreasing with both increasing 
y and increasing f t . The corrected values Q 2 COT and x cor were calculated in terms of ft cor 
and y cor using equation 

The correction was tested using the NC control sample (cf. section 3.3) by evaluating the 
means and widths of the distributions of S A = (A COI — A DA ) / A D a for A = Q 2 , x or y in all 
intervals of these variables (cf. section 4.2). Typically, the mean values were below 10% 
for x and Q 2 and below 5% for y; the r.m.s. widths were roughly 20% for x, 25% for Q 2 , 
and 10% for y. In the highest x and Q 2 bins, these estimates are based on only a few NC 
events and hence are of limited statistical significance. 

4.2 Choice of Intervals 

The kinematic ranges of the variables x, y and Q 2 covered by this analysis are x>0.006, 
0<y<l and Q 2 >200GeV 2 . For the measurement of differential cross sections in x and 
Q 2 , equal bins in log(x) and \og(Q 2 ) respectively were used (see tables ^, ^| and |]). The 
y range was subdivided in 5 equal intervals. These choices of bins provided roughly equal 
numbers of events in all intervals and ensured that the intervals were much wider than 
the experimental resolutions and the systematic shifts described in section 4.1. 



4.3 Background Subtraction 



The contamination from NC events, estimated using a lepto Monte Carlo sample |fl~9 



was negligible. The background from photoproduction processes was determined by ap- 
plying the selection algorithm to photoproduction events, simulated with herwig, which 
had transverse energy E t > 20 GeV (E t being determined from the momenta of all fi- 
nal state particles). While the simulated photoproduction events which passed the CC 
selection all had E t > 30 GeV, they were concentrated at low f. The estimated photo- 
production contamination in the lowest bins of x and Q 2 was ~ 15%. 
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No event passing the CC selection was found in the MC samples simulating heavy quark 
photoproduction, yielding an upper limit on the cross section for this background of less 
than 0.2% of the CC cross section. The heavy quark background is therefore neglected. 
No subtraction has been applied for non-ep background (cosmic and beam halo muons, p- 
gas reactions, and overlays of such events with ep reactions). This background is negligible 
after event selection and visual scanning (cf. section 3.2). 



4.4 Acceptance Correction and Unfolding 

Bin-to-bin migrations and the efficiencies of trigger and selection cuts were taken into ac- 
count using a bin-by-bin correction method, where the correction factors were calculated 
from the Monte Carlo simulation described in section 2. Typical acceptance corrections 
range from 40% to 80%, with an average of about 60% (cf. tables [l] and || to ||). The 
statistical errors of the correction factors were about 10 times smaller than the statistical 
errors of the data and have been neglected. In addition to the acceptance corrections, £, 
defined as the ratio of the number of events reconstructed in a given bin to the number 
of events generated in the same bin, tables || to ^| also show the purities, £, defined as the 
fraction of the events reconstructed in a given bin which were also generated in the same 
bin. 

The efficiency calculation was checked using the NC control sample. Both the overall 
acceptances and those of individual cuts (cf. table [I]) typically agreed with the MC values 
within 5%. 

The bin-by-bin unfolding procedure described above was cross checked by applying a 



single step matrix unfolding based on Bayes' theorem |3"3fl . The off-diagonal elements of 
the unfolding matrix were usually less than 25% of the corresponding diagonal elements, 
indicating that migrations between bins were small. In cases where the data are poorly 
modeled by the MC or where migrations play a dominant role, the two methods are ex- 
pected to yield significantly different results. The resulting differences were however found 
to be well below the statistical errors. They are included in the systematic uncertainties 
(cf. tables | and D. 



4.5 Systematic Uncertainties 

In the following, we summarize the studies which were performed in order to estimate the 
systematic uncertainties. For the case of da/dQ 2 , the dominant systematic uncertainties 
are detailed in tables |] and The remaining systematic errors were all below 5% and 
typically below 1%. Where appropriate, independent checks of systematic uncertainties 
have been calculated using the NC control sample. The systematic deviations observed 
in this sample agreed well with the values quoted below. 

• In the detector simulation, the calorimeter energy scale was varied by ±3% (±5% 
in the FCAL), corresponding to the level at which it is currently understood. Using 
the BAC it was checked that this variation also covers possible effects due to energy 
leakage out of the CAL. Using the BAC energy deposits to correct jet energies in 
CC events increased f t by about 2% on average. 

Such an energy scale uncertainty gives rise to uncertainties in the cross section which 
are typically between 5% and 10%, but can be up to 30% in the highest Q 2 and y 
bins. 



The systematic uncertainty of the photoproduction background subtraction was 
estimated as the quadratic sum of its statistical error and an additional ±50% 
uncertainty which accounted for differences between E t and ftt/E t spectra of pho- 
toproduction data and MC. 

The uncertainty of the background estimate mainly affected the lowest Q 2 and x 
bins. It was the dominant systematic error in the lowest x bin for e + p reactions. 

The reconstruction of the kinematic variables was modified in the following way: in 
each bin! of A = ft t or y the reconstructed values of f and y were multiplied with 
1/(1 + (5A)) (cf. section 4.1), where 5 A was averaged over the weighted NC control 
sample events in this bin. This corresponded to correcting on average to the DA 
values as derived from the NC control sample. Q 2 , x and y were calculated from 
the resulting values using equations [5| This changed the measured cross sections 
by 10% or less, except in the lowest and highest Q 2 bins, where variations of up to 
15% were seen. 

In order to check the simulation of the tracking, the combined acceptance of all cuts 
using tracking information was evaluated for the NC control sample, in each interval 
of Q 2 , x and y, both for data and for simulated NC events. The data/MC ratio 
of these acceptances was found to be consistent with unity. In almost all bins the 
deviation was less than 3%, which was used as the overall systematic uncertainty. 

A one step Bayes unfolding was used instead of the bin-by-bin unfolding (cf. dis- 
cussion in section 4.4). This variation caused uncertainties of 10 — 15% in the lowest 
and the highest bins of Q 2 , x and y (—25% in the highest y bin for e + p). 

In the trigger simulation, the first level trigger thresholds were increased to a level 
where the measured trigger efficiencies were close to 100%. Except in the lowest Q 2 
bins, where effects up to 4% were observed, this had a negligible effect indicating 
that the efficiency calculations were insensitive to the detailed shape of the trigger 
turn-on curves. 

In order to test the sensitivity to the details of the parton shower simulation, the 
meps model was used instead of the ARIADNE MC for calculating acceptance cor- 
rections and for unfolding. This affected the cross section results only negligibly. 

In addition, the selection cuts on ft, Acf), the track requirements, and the cuts at 
low ft t have been varied using the meps MC sample instead of the data, and the 
ariadne sets for acceptance corrections and unfolding. This tested whether the 
simulated shapes of the relevant distributions differ in the vicinities of these cuts. 

Each of the cut variations changed the results by less than typically 1%, except the 
ft variation in the lowest x and Q 2 bins (up to 5%) and the track requirement (up 
to 3% in some bins). 

Different parton distribution functions (MRSD0 ||, MRSD'_ gjj, GRV 11 and 



CTEQ 2pM |J7]]) were used in the MC samples used for the acceptance corrections.0 
The resulting variations were small and have been neglected. 



9 For f t , the bin boundaries were 0-20-30-40-50-oo GeV. 

10 The most significant variation of the Standard Model cross section prediction due to parton distri- 
butions is seen in da/dx at 0.01 ^ x 0.1 (cf. section 5.1). 
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The overall systematic uncertainty for each result (including ratios of cross sections and 
the W mass) was calculated by evaluating the variation of this result due to each of 
the modifications of the analysis procedure described above (including the luminosity 
uncertainties of ±2.0% for the e + p sample and ±3.0% for the e~p sample) and then 
separately summing in quadrature the positive and negative deviations. 

5 Results 

The results will be presented in three sections: cross sections and the determination 
of Myy, jet rates and profiles, and finally the search for events with a large rapidity 
gap. The abbreviation 'SM' (for Standard Model) is used to indicate the theoretical 
predictions. Unless otherwise noted, these were estimated using lepto and the MRSA 
parton distributions. The propagator four-momentum was used to calculate kinematic 
quantitites. 

5.1 Cross Sections and Cross Section Ratios 

Differential cross sections as functions of Q 2 , x and y 

The values of the cross section, <jj, in bins of Q 2 , for e + p and e~p data, are listed in 
tables H and respectively. Also shown are the dominant systematic errors as described 
in section 4.5. The results are plotted in figure |2|a, where the SM predictions are also 
shown. 

The results for Q 2 >200 GeV 2 in bins of x and y are given in tables |] and |5[ The differential 
CC cross sections are shown in figures 0a and £|a together with the corresponding SM 
predictions. 

Figures 0b, |3|b and f|b show the ratios <r e + p ->vx I ff e -p->vx as functions of the respective 
variables. The curves again represent the SM predictions. 
Several observations can be made: 

• No significant deviations of the measured cross sections from the SM predictions are 
observed. This is quantified by the \ 2 values which are listed in table ||. 

• The differential cross sections da/dQ 2 fall steeply with Q 2 (cf. figure 0), reflecting 
the influence of the W propagator, the decrease of the parton densities with in- 
creasing x = Q 2 /sy, and the (1 — y) 2 terms in the cross section. The differences 
between the e + p and e~p cross sections are due to the different x and y behavior 
(cf. discussion below). 

• As can be seen from figure ^|, at low x, da/dx is about the same for e + p and e~p 
scattering. This reflects the fact that the cross sections (equations |l]|J) become equal 
if valence quark distributions can be neglected and if d(x, Q 2 ) = u(x, Q 2 ). 

The decrease of da/dx with increasing x is more rapid for e + p scattering than for 
e~p scattering. This behavior is also expected from equations [l] and ||: in e + p 
reactions, scattering on valence quarks is reduced by the (1 — y) 2 factor; in addition, 
the e + couples to d(x, Q 2 ), which is suppressed at high x relative to u(x, Q 2 ) which 
is relevant for e~p scattering (cf. e.g. [|38|j). 
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The x dependence of the e + p to the e p cross section ratio is illustrated by the 
ratios for x < 0.1 and for x > 0.1 (with Q 2 > 200 GeV 2 ): 



r o.9i 



+0.41 
-0.27 



[stat] ± 0.04 [syst] 



for < x < 0.1 



= < 



0.26 



+0.14 
-0.09 



[stat] ± 0.01 [syst] 



for 0.1 < x < 1 . 



The corresponding SM predictions are 0.64 and 0.32. The measured ratios confirm 
the suppression of the e + p CC cross section at high x. 

The SM predictions for da/dx have also been evaluated using different parton distri- 
bution functions (MRSD'_, GRV, CTEQ 2pM). The largest differences have been ob- 
served between the sets MRSD'_ and CTEQ 2pM for x between 0.01 and 0.1, where 
the cross section predicted by MRSD'_ is up to 10% lower than the CTEQ 2pM val- 
ues, both for e + p and e~p scattering. The MRSA, CTEQ 2pM and GRV predictions 
agree to better than 5%. 

• The measured differential cross sections da/dy agree with the theoretical predic- 
tions. The SM curves in figure |] reveal that the expected shapes are very similar 
for e + p and e~p scattering, indicating that the relative contributions of the (1 — y) 2 
terms are about equal in both cases. 

As y — > 1, the cross section ratio e + p/e~p is given by 1/ ( - - + lj (assuming 
u = u s and c = c, the indices v and s indicating valence and sea quark distributions, 
respectively). The fact that the prediction, corresponding to u v /{u s + c) « 1, is 
above the data for y > 0.4, is only moderately significant in view of the large errors. 

At y — > 0, the cross section ratio e + p/e~p approaches (d v + S)/(u v + S), where 
S = d s +u s + s + c (assuming again that the quark and antiquark sea distributions 
are identical for each flavor). 

In order to further study the influence of the V—A helicity structure and the parton 
distributions, figure || shows da/dy separately for x < 0.1 (where sea quarks dominate 
the cross section) and for x > 0.1 (where the valence contribution is larger). Here three 
equal size y bins are used. 

For low x, the predicted y shapes are similar for e + p and e~p scattering, the differences 
being mainly due to residual valence contributions. The most significant deviation of the 
measured cross sections from the prediction is observed in the lowest y bin (y < 0.33), 
where the measured a e - p ^ v x is about two standard deviations below the SM prediction. 

For x > 0.1, both the e + p and the e~p cross sections decrease with increasing y. For e~p, 
this can be attributed to residual sea quark contributions and to the W propagator which 
enters via the dependence of the average Q 2 on y. For e + p scattering, the decrease of 
da/dy with increasing y is even faster due to the (1 — y) 2 term in the cross section. All 
measured values are compatible with the predictions, indicating that the slight deviations 
observed at y ~ 0.3 in the e + p/e~p ratios are due to the low-x data. 
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Integrated cross sections for Q 2 > 200 GeV 2 



The integrated cross sections for Q 2 > 200 GeV 2 are obtained by summing the cross 
sections in the Q 2 bins of tables |2| and ^ resulting in§ 

a e+p ^ vx {Q 2 > 200 GeV 2 ) = 30.3 ± \\ [stat] ±™ [syst] pb 
a e ~ p ^ uX {Q 2 > 200 GeV 2 ) = 54.7 t 1 ™ [stat] ±H [syst] pb . 

These measurements are compatible with the respective SM predictions of 32.3 pb and 
65.8 pb. The e + p/e~p ratio is 

^^(Q 2 > 200 GeV 2 ) = 0.55 ±£2 [stat] i°;° 2 [syst] , (6) 

in agreement with the SM prediction of 0.49. 
Determination of the W mass 

The Q 2 dependence of the CC cross section is largely determined by the propagator term 
V(M W , Q 2 ) = [M W /(M W + Q 2 )} 2 . In order to determine the value of M w , the differential 
cross section was factorized according to dcr/dQ 2 = V(Mw,Q 2 ) • &(Q 2 )- The function 
$(<5 2 ), containing the Q 2 dependence of the parton densities and couplings, was taken 
from the MC simulation, using the MRSA parton distributions. A binned log-likelihood 
fit with M\y as the free parameter, applied simultaneously to the e + p and e~p data 
samples, yielded 

M w = 79 t 8 7 [stat] t\ [syst] GeV , (7) 



in agreement with the average value Mw = 80.22 ± 0.26 GeV [5{| obtained from direct 
measurements of the W mass at pp colliders || and also with the recent HI measurement 
. The systematic uncertainty of M w has been evaluated using the method described in 
section 4.5. The dominant contribution is the calorimeter energy scale. 



5.2 Jet Analysis 

An (77, 0)-cone jet finding algorithm [HD| has been applied to the data, using a cone 



radius, R = (A0 2 + Ar/ 2 )^ of 0.7. Pre-clusters are formed around calorimeter cells 
with transverse energies larger than 0.3 GeV, and the final clusters are called jets if their 
transverse energies exceed 6 GeV and their pseudorapidities 7] ]Ct are less than 2.5 (i.e. polar 
angles greater than 9°). This jet analysis is done in the laboratory system. A detailed 
analysis of jets in NC reactions, covering the region Q 2 ^ 10 3 GeV 2 , has previously been 
published in [p| . 

The results of the present jet analysis turn out to be insensitive to the choice of the cone 
radius (R — 1 and R = 0.5 have also been tested). The constraint on rj set was applied 
in order to avoid effects from the fragmentation of the proton remnant. In fact, the 
proton remnant manifests itself as a moderate enhancement of the transverse energy flow 
around the forward edge of acceptance, and for 7]i et = 2.5 this enhancement is still clearly 
separated from the E t flow associated with the jet. 



11 Summing the single cross sections implies using an average of the acceptances evaluated for the single 
bins, which is slightly different from the global acceptance quoted in table 0. 
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For all jets found, the distributions of the jet momentum transverse to the beam axis, pj e , 
is shown in figures ||a and |7|a for e + p and e~p data, respectively, pf has been corrected 
for energy loss in the inactive material of the detector^ by using the average difference 
of measured and true values of pl et as determined in the MC simulation as a function of 
the measured p f jet . Note that the pl et spectrum is much harder for the e~p than for the 
e + p data, as is expected from the harder Q 2 distributions. Figures |6|b and [7|b show the 
distributions of r]^ et . The transverse energy flows measured in pseudorapidity (At/) and 
azimuthal angle (A</>) relative to the jet axis, are plotted in figures |6|c,d and |7|c, d. The 
hadronic energy flow produced between the jet and the proton remnant is clearly observed 
in the excess of transverse energy at high An, i.e. towards the proton remnant. 
All these figures also show the predictions of the ARIADNE MC. The pl et and rj ict dis- 
tributions are also compared to the meps simulation. The MC distributions have been 
normalized to the number of jets observed in the data. For the E t flow distributions, 
the differences between ARIADNE and meps (not shown) are slight. Good agreement is 
observed between the data and both the ARIADNE and the meps predictions, indicating 
in particular that the jet properties are well described by the MC models. 
The jet multiplicity distributions are summarized in table [7[ All selected events have at 
least one jet. Events with two or more jets are expected from hard QCD subprocesses, in 
particular from boson-gluon fusion (BGF). The MC predictions for the jet multiplicities, 
which include the 0(a s ) matrix elements, exhibit significant differences, amounting to 
nearly a factor of two in the 2-jet rates. However, within their statistical errors, the data 
are compatible with both models, favoring values intermediate between ARIADNE and 

MEPS. 



5.3 Search for Large Rapidity Gap Events 

Events with a large rapidity gap (LRG) between the outgoing proton system and the 
rest of the hadronic final state have been observed at HERA in neutral-current DIS f7|. 
These events are generally understood to be of a diffractive nature and to result from the 
exchange of a colorless object (usually called the pomeron) with the quantum numbers 
of the vacuum. Evidence for a partonic structure of the pomeron in DIS has also been 
observed by the HERA experiments ||. 

In CC processes, the coupling of the exchanged W is sensitive to the flavor of the pomeron 
constituents, which could provide additional information on the pomeron structure. The 
search also is sensitive to high Q 2 diffractive production of exclusive hadronic states such 
as vector or axial-vector mesons, complementing corresponding NC studies . 



Similar to the prescription used in a previous NC analysis 0, LRG events are identified 
using two variables. ?7 max and Oh- f] maiX is defined as the maximum pseudorapidity of any 
calorimeter clusterta or cell with energy greater than 400 MeV. The global quantity Oh is 
given by the energy weighted mean polar angle of the energy deposits in the calorimeter, 
cos 9h = J2i P l z / J2i E 1, , where the sums run over all calorimeter cells. In the naive quark 
parton model, Oh is the scattering angle of the struck quark. 

Figure || shows a scatter plot of cos Oh versus r/ max for the combined e + p and e~p samples. 
The rectangle indicates the region ?? max <2.5 and cos#j/<0.75, which is taken from M as 



12 This was the only correction applied in this jet analysis. 

13 Calorimeter clusters are groups of adjacent cells which have energy deposits above the noise 
thresholds. 
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the definition of LRG events ("LRG requirement"). One event of the e + p sample, shown 
in figure passes these LRG cuts. The reconstructed kinematic variables for this event 
are #=14 ± 2 GeV, Q 2 =300 ± 70 GeV 2 , x=0.0093 ± 0.0015 and y=0.35 ± 0.10. 

In order to estimate event selection efficiencies and expected event rates in the e + p case, 
the LRG analysis was also applied to the following CC MC sets (at the generator level, 
^max was defined as the maximum pseudorapidity of all particles with momenta greater 
than 400 MeV and with rj < 4.5, where the latter requirement excluded the proton in 
diffractive reactions): 

• ariadne. This set uses lepto for the differential cross section and includes the 
simulation of diffractive-like processes via soft color interactions between the per- 
turbatively produced partons and color-charges in the proton remnant E31. The 



fraction of events with r] max < 2.5 at the generator level was found to be 1.3%. 

• meps. This set also is based on lepto, but does not simulate soft color interactions. 
Here the rapidity gap between proton direction and hadronic system is exponentially 
suppressed. At the generator level, 0.45% of the events had rj max < 2.5. Most of 
the events passing the LRG requirement had a smaller rapidity gap at the generator 
level. 

• pompyt. This generator simulates only diffractive reactions and is based on a 
factorizable model for high-energy diffractive processes where, within the pythia 
framework, the proton emits a pomeron, whose constituents take part in a hard 
scattering process with the W boson. For this analysis a hard quark density distri- 
bution in the pomeron was assumed (oc (3(1-/3), where j3 is the momentum fraction 
of the quark relative to the pomeron momentum), which in addition satisfied the 
momentum sum rule.0 

It is interesting to note that, according to the pompyt simulation, only 18% of all 
CC diffractive events with Q 2 > 200 GeV 2 have r] mauX < 2.5 at the generator level, 
41% of which pass the CC selection and fulfill the LRG requirement. 

The fractions of LRG events are summarized in table § (for the e + p NC control sample this 
ratio was determined using the event weights as described in section 3.3). They roughly 
agree between the CC data, the ariadne simulation and the NC control sample, whereas 
for the meps sample, the fraction of LRG events is below that predicted in ariadne. 
The selection efficiency £lrg> defined as 

no. of events with Q 2 or > 200 GeV 2 which pass CC cuts and LRG requirement 
Lh( ' no. of events which have Q 2 > 200 GeV 2 and ^ max < 2.5 at generator level ' 

is also shown in table [5] for the three MC sets. Using the ariadne value of £lrg = 0.43, 
the cross section for e + p CC interactions with ?7 max < 2.5 was calculated to be 

Ve+ P -,vx(Q 2 > 200 GeV 2 ; Vmax < 2.5) = 0.8 [stat] ± 0.1 [syst] pb. (8) 

The systematic uncertainty of this result was derived from the difference of the efficiencies 
£lrg obtained using the meps and pompyt sets. 



14 The latter assumption conflicts with measurements || of the NC diffractive cross section. Hence 
pompyt predictions of absolute cross sections were not used. 
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6 Summary and Conclusions 



The cross sections for the deep inelastic charged-current interactions e + p — > VX and 
e~p — > vX have been measured in the kinematic region x > 0.006 and Q 2 > 200 GeV 2 , 
using the e + p and e~p data collected with the ZEUS detector during 1993 and 1994. For 
the integrated cross sections, 



are obtained, in good agreement with the Standard Model expectations of 32.3 pb and 
65.8 pb. 

The differential cross sections have been measured as functions of Q 2 , x and y. The 
differences between e + p and e~p scattering, which are expected from the V— A helicity 
structure of the weak interaction and from the quark content of the proton, are clearly 
observed. No significant deviation from the Standard Model predictions is found. 

The measured Q 2 dependence of the differential cross sections for e + and e~ induced 
CC D1S reactions is used to determine the W mass. The measured value of My/ is 
79 If 1 4 GeV, in good agreement with direct measurements performed at pp colliders. 

A jet analysis has been performed and CC events with multiple jets have been observed. 
The distributions of the jet transverse momentum and rapidity as well as transverse energy 
flow jet profiles have been measured and are compared to the 0(a s ) predictions of the 
ARIADNE and meps Monte Carlo models. Good agreement is observed in both the jet 
profiles and the jet rates. 

A search for charged-current events with a large rapidity gap between the observed 
hadronic system and the outgoing proton system yielded one candidate event in the 
e + p sample, corresponding to a cross section of a e + p ^x(Q 2 > 200 GeV 2 ; i] max < 2.5) = 
0.8 ±5:f [stat] ±0.1 [syst]pb. 



Acknowledgments 

We appreciate the contributions to the construction and maintenance of the ZEUS de- 
tector of the many people who are not listed as authors. The HERA machine group and 
the DESY computing staffs are especially acknowledged for their efforts to provide excel- 
lent operation of the collider and the data analysis environment. We thank the DESY 
directorate for strong support and encouragement. 



a e+p _ tVX (Q 2 >200GeV 2 ) 
a e - p ^ uX {Q 2 >200GeV 2 ) 



30.3 1 H [stat] t\l [syst] pb 
54.7™ [stat] ±§;| [syst] pb 



15 



References 

[1] CDHS Collaboration, H.Abramowicz et al., Z. Phys. C25(1984)29; 
CDHSW Collaboration, J.P.Berge et al., Z. Phys. C49(1991)187; 
CCFR Collaboration, E.Oltman et al, Z. Phys. C53(1992)51; 
BEBC WA21 Collaboration, G.T.Jones et al, Z. Phys. C62(1994)575. 

[2] HI Collaboration, T.Ahmed et al., Phys. Lett. B324(1994)241. 

[3] ZEUS Collaboration, M.Derrick et al, Phys. Rev. Lett. 75(1995)1006. 

[4] HI Collaboration, S.Aid et al., Z. Phys. C67(1995)565. 

[5] HI Collaboration, S.Aid et al., DESY 96-046 (1996). 

[6] CDF Collaboration, F.Abe et al., Phys. Rev. Lett. 75(1995)11; 
CDF Collaboration, F.Abe et al., Phys. Rev. D52(1995)4784; 

DO Collaboration, S.Abachi et al., Proc. Topical Workshop on Proton-Antiproton 

Collider Physics, Tsukuba, Japan (1993); 

UA2 Collaboration, J.Alitti et al., Phys. Lett B241(1990)150. 

[7] ZEUS Collaboration, M.Derrick et al., Phys. Lett. B315(1993)481; 
ZEUS Collaboration, M.Derrick et al., Phys. Lett. B332(1994)228; 
HI Collaboration, T.Ahmed et al., Nucl. Phys. B429(1994)477. 

[8] ZEUS Collaboration, M.Derrick et al., Phys. Lett. B346(1995)399; 
HI Collaboration, T.Ahmed et al., Phys. Lett. B348(1995)681; 
ZEUS Collaboration, M.Derrick et al., Phys. Lett. B356(1995)129. 

[9] ZEUS Collaboration, M.Derrick et al., Z. Phys. C68(1995)569. 

[10] ZEUS Collaboration, M.Derrick et al., DESY 96-018 (1996). 

[11] ZEUS Collaboration, M.Derrick et al., Z. Phys. C65(1995)379. 

[12] ZEUS Collaboration, M.Derrick et al., Phys. Lett. B293(1992)465. 

[13] The ZEUS Detector, Status Report 1993, DESY (1993). 

[14] M.Derrick et al., Nucl. Inst. Meth. A309(1991)77; 
A.Andresen et al., Nucl. Inst. Meth. A309(1991)101; 
A.Bernstein et al., Nucl. Inst. Meth. A336(1993)23; 

A. Caldwell et al., Nucl. Inst. Meth. A321(1992)356. 

[15] H.Abramowicz et al., Nucl. Inst. Meth. A313(1992)126. 
[16] C.Alvisi et al., Nucl. Inst. Meth. A305(1991)30. 

[17] N.Harnew et al., Nucl. Inst. Meth. A279(1989)290; 

B. Foster et al., Nucl. Phys., Proc. Suppl. B32(1993); 
B.Foster et al., Nucl. Inst. Meth. A338(1994)254. 

[18] J.Andruszkow, DESY 92-066 (1992). 



16 



[19] lepto 6.3: G.Ingelman, Proc. 1991 Workshop on Physics at HERA, ed. 
W.Buchmiiller and G.Ingelman (DESY, Hamburg, 1992), Vol. 3, p. 1366. 

[20] Heracles 4.4: A.Kwiatkowski, H.Spiesberger and H.J.Mohring, Proc. 1991 Work- 
shop on Physics at HERA, ed. W.Buchmiiller and G.Ingelman (DESY, Hamburg, 
1992), Vol. 3, p.1294. 

[21] django 6.1: G.Schuler and H.Spiesberger, Proc. 1991 Workshop on Physics at 
HERA, ed. W.Buchmiiller and G.Ingelman (DESY, Hamburg, 1992), Vol. 3, p. 1419. 

[22] A.D.Martin, R.G.Roberts, W.J.Stirling, Phys. Rev. D50(1994)6734. 

[23] ariadne 4.06: L.Lonnblad, LU TP-89-10; 
L.Lonnblad, Comp. Phys. Comm. 71(1992)15. 

[24] pythia 5.7 and jetset 7.4: T.Sjostrand, CERN-TH 7112-93 (1994); 
T.Sjdstrand, LU-TP-95-20 (1995). 

[25] herwig 5.8: B.R.Webber, Proc. 1991 Workshop on Physics at HERA, ed. 
W.Buchmiiller and G.Ingelman (DESY, Hamburg, 1992), Vol. 3, p. 1363. 

[26] aroma 2.1: G.Ingelman and G.Schuler, Proc. 1991 Workshop on Physics at HERA, 
ed. W.Buchmiiller and G.Ingelman (DESY, Hamburg, 1992), Vol. 3, p.1346. 

[27] pompyt 1.0: P.Bruni and G.Ingelman, Proc. Europhysics Conference on HEP, Mar- 
seilles, France (1993), p.595. 

[28] G.Ingelman and PSchlein, Phys. Lett. B152(1985)256. 

[29] R.Brun et al., CERN DD/EE-84-1 (1987). 

[30] F.Jacquet and A.Blondel, Proceedings of the study for an ep facility for Europe, 
DESY 79-48 (1979), p.391. 

[31] H.Abramowicz, A.Caldwell, R.Sinkus, Nucl. Inst. Meth. A365(1995)508. 

[32] S.Bentvelsen, J.Engelen, P.Kooijman, Proc. 1991 Workshop on Physics at HERA, 
ed. W.Buchmiiller and G.Ingelman (DESY, Hamburg, 1992), Vol. 1, p.23. 

[33] G.D'Agostini, DESY 94-099 (1994). 

[34] A.D.Martin, R.G.Roberts, W.J.Stirling, Phys. Rev. D47(1993)867. 

[35] A.D.Martin, R.G.Roberts, W.J.Stirling, Phys. Lett. 306B(1993)147; 
erratum in Phys. Lett. 309B(1993)492. 

[36] M.Gluck, E.Reya, A.Vogt, Z. Phys. C53(1992)127. 

[37] CTEQ Collaboration, J.Botts et al, Phys. Lett. B304(1993)159. 

[38] BEBC WA21 Collaboration, G.T.Jones et al., Z. Phys. C62(1994)601. 

[39] L.Montanet et al, Phys. Rev. D50(1994)1173. 



17 



[40] UA1 Collaboration, G.Arnison et al, Phys. Lett. B123(1983)115; 

J.Huth et al., Proc. 1990 DPF Summer Study on High-Energy Physics, Snowmass, 
Colorado, ed. E.L.Berger (World Scientific, Singapore, 1990, p. 134). 

[41] ZEUS Collaboration, M.Derrick et al, Z. Phys. C67(1995)81. 

[42] ZEUS Collaboration, M.Derrick et al, Phys. Lett. B356(1995)601; 
HI Collaboration, S.Aid et al., DESY 96-023. 

[43] A.Edin, G.Ingelman, J.Rathsman, DESY 95-145 (1995). 



18 



cut 
description 


e + p 


e p 


data 


MC(Q 2 >200 GeV 2 ) 


data 


MC(Q 2 >200 GeV 2 ) 


Q 2 > 200 GeV 2 (MC only) 




100 % 




100 % 


Trigger and preselection 


6307 


83.2% 


1450 


86.1% 


ft > 11 GeV 


2981 


80.0% 


940 


83.1% 


Z vtx < 45 cm 


1283 


71.0% 


500 


74.7% 


rry ft > o.6 


958 


70.3% 


301 


73.9% 


good track 


741 


68.2% 


237 


71.7% 


|A0| < 1 


342 


66.8% 


113 


70.3% 


vjb < 0.8 


279 


66.2% 


98 


69.2% 


Ft/E t > 0.4 


208 


65.0% 


59 


68.1% 


If***/ ft > 0.1 


87 


63.0% 


36 


66.3% 


Muon, timing, sparks 


59 


63.0% 


30 


66.3% 


Visual scan (data only) 


56 




30 





Table 1: Data reduction by the selection cuts described in section 3.2. The columns 
denoted "data" indicate the numbers of events remaining after applying the cuts. The 
"MC" columns show the cumulative selection efficiencies for the ARIADNE MC simulation. 
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4- 
e~p 


Q 2 [GeV 2 ] 


200 
437 


437 
955 


955 
2089 


2089 
4517 


4517 
10000 




N cvt 


8 


15 


12 


14 


7 


56 


(Q 2 ) [GeV 2 ] 


oil) 


(v 


14oU 


oil ([) 


oolU 


1 m n 

lyiu 


AT 


1 9 


n 9 


n 

u 


n 

u 


n 

u 


1 A 


t 
s 


0.38 


0.58 


0.70 


0.76 


0.71 


0.62 




0.76 


0.72 


0.73 


0.76 


0.79 


1.0 


<j exp [pb] 


6.1 


8.7 


5.8 


6.3 


3.4 


30.3 


A<j cxp [stat] [pb] 


+ 3.0 
- 2.2 


+ 2.9 
— 2.3 


+ 2.2 
- 1.7 


+ 2.2 
- 1.7 


+ 1.8 
- 1.3 


+ 5.5 
- 4.2 


A(t cx p [syst] [pb] 


+ 1.3 

- 2.7 


+ 0.6 
- 0.6 


+ 0.3 
- 0.3 


+ 0.6 
- 0.3 


+ 0.7 

- 0.4 


+ 1.6 
- 2.6 


(j th [pb] 


6.2 


8.3 


8.3 


6.4 


2.5 


32.3 


Systematic checks 


Check 


A(7 ex P [%] 


E scale +3% (+5% in FCAL) 


12.6 


6.3 


2.1 


- 3.3 


-10.8 


2.1 


E scale -3% (-5% in FCAL) 


-18.7 


- 3.0 


- 3.3 


3.8 


13.6 


- 2.0 


Unfolding 


- 0.9 


1.1 


1.7 


- 1.5 


13.8 


-- 1.0 


Q 2 ,x,y reconstruction 


-15.7 


- 2.1 


2.0 


7.6 


5.8 


- 0.3 


A bg + AA bg 


-36.2 


- 3.4 


0.0 


0.0 


0.0 


- 7.3 


A bg - AA bg 


16.1 


1.5 


0.0 


0.0 


0.0 


3.3 



Table 2: Cross sections for CC e + p scattering in intervals of Q 2 , and for Q 2 > 200 GeV 2 . 
The rows denoted A evt , (Q 2 ), A bg , (, £, <r exp , and a th show the raw event numbers, the 
mean Q 2 according to the SM prediction, the estimated number of background events from 
photoproduction, the bin-by-bin acceptance corrections, the purities, the experimental 
cross sections, and the theoretical cross sections, respectively. For <r cxp the statistical 
errors Acr exp [stat] and systematic uncertainties Acr exp [syst] are quoted. The latter are 
calculated from the deviations resulting from the systematic studies described in the text 
(including the luminosity uncertainty of 3% and the overall uncertainty of 3% assigned 
to the efficiency of those selection cuts involving track quantities). The dominant contri- 
butions are detailed in the lower part of the table. 
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e p 


Q 2 [GeV 2 ] 


200 
437 


437 
955 


955 
2089 


2089 
4517 


4517 
10000 


10000 
21877 


> 200 


N evt 


3 


3 


7 


6 


9 


2 


30 


(Q 2 ) [GeV 2 ] 


310 


670 


1470 


3140 


6700 


13900 


3340 


N hg 


0.3 


0.1 














0.4 


C 


0.39 


0.56 


0.74 


0.78 


0.74 


0.59 


0.66 




0.74 


0.73 


0.70 


0.74 


0.80 


0.83 


1.0 


(T cxp [pb] 


8.4 


6.3 


11.6 


9.4 


14.9 


4.1 


54.7 


A(T cxp [stat] [pb] 


+ 8.2 
- 4.6 


+ 6.2 
- 3.4 


+ 6.3 
- 4.3 


+ 5.6 
- 3.7 


+ 6.8 
- 4.9 


+ 5.4 
- 2.7 


+15.9 
- 9.8 


A(j cx p [syst] [pb] 


+ 1.6 

- 2.7 


+ 0.6 
- 0.7 


+ 0.7 
- 0.7 


+ 0.8 
- 0.6 


+ 1.9 

- 1.6 


+ 1.5 
- 1.2 


+ 2.8 
- 2.9 


a th [pb] 


7.9 


12.2 


15.1 


15.7 


10.6 


3.7 


65.8 


Systematic checks 


Check 


Aa cxp [%] 


E scale +3% (+5% in FCAL) 


14.2 


7.1 


2.8 


0.9 


- 3.9 


-24.9 


2.0 


E scale -3% (-5% in FCAL) 


-10.9 


- 7.8 


- 2.6 


1.2 


9.5 


26.9 


1.1 


Unfolding 


- 0.5 


0.0 


- 3.3 


- 4.2 


- 9.3 


24.3 


- 2.1 


Q 2 ,x,y reconstruction 


-15.2 


- 5.9 


2.2 


6.8 


7.8 


-14.5 


- 1.0 


N hg + AN hg 


-25.9 


- 5.0 


0.0 


0.0 


0.0 


0.0 


- 4.0 


N hg - AN hg 


11.1 


2.7 


0.0 


0.0 


0.0 


0.0 


1.7 



Table 3: Cross sections for CC e~p scattering in intervals of Q 2 , and for Q 2 > 200 GeV 2 . 
The rows denoted N evt , (Q 2 ), N hg , (, £, <r cxp , and a th show the raw event numbers, the 
mean Q 2 according to the SM prediction, the estimated number of background events from 
photoproduction, the bin-by-bin acceptance corrections, the purities, the experimental 
cross sections, and the theoretical cross sections, respectively. For <r cxp the statistical 
errors Acr exp [stat] and systematic uncertainties Acr exp [syst] are quoted. The latter are 
calculated from the deviations resulting from the systematic studies described in the text 
(including the luminosity uncertainty of 2% and the overall uncertainty of 3% assigned 
to the efficiency of those selection cuts involving track quantities). The dominant contri- 
butions are detailed in the lower part of the table. 
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1 

e~p 


X 


0.006 
0.014 


0.014 
0.033 


0.033 
0.076 


0.076 
0.17 


0.17 
0.40 


0.40 
1.0 


N cvt 


3 


20 


11 


13 


9 





(x) 


0.010 


0.023 


0.052 


0.115 


0.249 


0.480 


N hg 


1 /I 
1.4 


n 
u 


n 
u 


n 
u 


n 
u 


n 
u 


/- 
s, 


n a 1 

U.41 


U. / o 


u. < y 


U.Oo 


n a s 
U.4o 
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t 
K 


0.68 


0.72 


0.77 


0.82 


0.88 


0.84 


a cxp [pb] 


1.3 


8.7 


4.7 


6.5 


6.4 


0.0 


Aa exp [stat] [pb] 


+ 1.3 
- 0.7 


+ 2.4 
- 1.9 


+ 1.9 
- 1.4 


+ 2.3 
- 1.8 


+ 2.9 
- 2.1 


+ 3.5 
- 0.0 


Aa exp [syst] [pb] 


+ 1.0 

- 1.3 


+ 0.5 
- 0.7 


+ 0.4 
- 0.2 


+ 0.5 
- 0.7 


+ 0.8 
- 0.6 


+ 0.0 
- 0.0 


a th [pb] 


3.4 


6.3 


8.4 


8.3 


4.2 


0.5 




e p 


X 


0.006 
0.014 


0.014 
0.033 


0.033 
0.076 


0.076 
0.17 


0.17 
0.40 


0.40 
1.0 


iVcvt 


3 


5 


5 


9 


8 





(x) 


0.010 


0.024 


0.053 


0.119 


0.256 


0.489 




0.4 

















C 


0.41 


0.77 


0.84 


0.72 


0.55 


0.17 




0.61 


0.73 


0.73 


0.85 


0.87 


0.83 


a cxp [pb] 


7.8 


7.9 


7.3 


15.3 


17.6 


0.0 


Aa exp [stat] [pb] 


+ 7.6 
- 4.2 


+ 5.3 
- 3.4 


+ 4.9 
- 3.2 


+ 7.0 

- 5.1 


+ 8.7 
- 6.3 


+13.2 
- 0.0 


Aa cxp [syst] [pb] 


+ 1.5 
- 2.5 


+ 0.7 
- 0.7 


+ 0.4 
- 0.5 


+ 1.1 

- 2.0 


+ 3.0 
- 1.6 


+ 0.0 
- 0.0 


a th [pb] 


4.2 


8.4 


14.9 


19.2 


14.8 


2.7 



Table 4: Cross sections for CC e + p (top) and e~p (bottom) scattering for Q 2 > 200 GeV 2 
in intervals of x. The rows denoted A^ cvt , (x), N^, g , (, £, <r exp , and a th show the raw event 
numbers, the mean x according to the SM prediction, the estimated number of background 
events from photoproduction, the bin-by-bin acceptance corrections, the purities, the 
experimental cross sections, and the theoretical cross sections, respectively. For <r exp the 
statistical errors Aa cxp [stat] and systematic uncertainties A<r cxp [syst] are quoted. 
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0.79 


0.77 


0.70 


0.80 


(j cxp [pb] 


9.6 


10.7 


4.2 


1.4 


2.3 


Aa exp [stat] [pb] 


+ 3.3 
- 2.5 


+ 2.5 
- 2.1 


+ 1.7 

- 1.3 


+ 1.4 

- 0.8 


+ 5.4 
- 1.9 


Aa exp [syst] [pb] 


+ 0.5 
- 0.7 


+ 0.6 
- 0.5 


+ 0.6 
- 1.0 


+ 0.2 
- 0.1 


+ 0.6 
- 1.1 


a th [pb] 


11.8 


8.3 


5.4 


3.7 


3.0 






y 




0.0 


0.2 


0.4 


0.6 


0.8 


e p 












0.2 


0.4 


0.6 


0.8 


1.0 




6 


5 


11 


5 


3 


(y) 


0.101 


0.293 


0.496 


0.696 


0.895 







0.1 


0.3 








C 


0.51 


0.91 


0.81 


0.77 


0.20 




0.92 


0.79 


0.81 


0.74 


0.90 


(j cxp [pb] 


14.4 


6.6 


16.2 


7.9 


19.0 


Aa exp [stat] [pb] 


+ 8.6 
- 5.7 


+ 4.4 
- 2.9 


+ 6.4 
- 4.9 


+ 5.2 
- 3.4 


+18.5 
-10.3 


Aa exp [syst] [pb] 


+ 1.2 
- 1.6 


+ 0.3 
- 0.4 


+ 1.3 
- 1.4 


+ 0.8 
- 0.5 


+ 4.4 
- 7.6 


a th [pb] 


23.7 


15.6 


11.0 


8.7 


6.9 



Table 5: Cross sections for CC e + p (top) and e~p (bottom) scattering for Q 2 > 200 GeV 2 
in intervals of y. The rows denoted iV evt , (y), iVb g , C> £> <r exp , and <r th show the raw event 
numbers, the mean y according to the SM prediction, the estimated number of background 
events from photoproduction, the bin-by-bin acceptance corrections, the purities, the 
experimental cross sections, and the theoretical cross sections, respectively. For <r exp the 
statistical errors A<r cxp [stat] and systematic uncertainties A<r cxp [syst] are quoted. 
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Distribution 


x 2 


-^DOF 


P(X 2 , N DOF ) 


da/dQ 2 (e + p) 


1.9 


5 


0.86 


da/dQ 2 (e'p) 


3.2 


6 


0.78 


da/dx (e + p) 


9.9 


5 


0.08 


da/dx (e~p) 


3.7 


5 


0.59 


da/dy (e + p) 


5.0 


5 


0.42 


da/dy (e~p) 


7.8 


5 


0.17 



Table 6: Quality of agreement between experimental results and SM predictions: for the 
various differential cross sections \ 2 ■, ^dof, and p(x 2 , -^dof) denote the x 2 value calcu- 
lated using statistical errors only, the number of degrees of freedom, and the corresponding 
^-probability. 



24 
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A'jets 
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ARIA. 
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A^evt 
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pMC 

"""jet 
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0.4 ±0.1 


0.4 ±0.1 





o.ol8:8 


0.3 ±0.1 


0.4 ±0.1 
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87.5 til 


90.9 ±0.4 
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91.6 ±0.4 
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8.4 ±0.4 
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3 
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7.8 ±0.4 
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1 


l-8±tg 


0.3 ±0.1 


1.1 ±0.2 





u - u -0.0 


0.3 ±0.1 


1.2 ±0.2 



Table 7: Jet rates for e+p and e~p data. The columns denoted iV ev t show the numbers of 
events in the CC samples which have iVj ets jets in addition to the proton remnant (for the 
jet definition see text). The Rf^ a are the corresponding jet rates, which are also given 
for the ARIADNE and MEPS MC simulations (i?H c ). Only statistical errors are shown. 



data set 


-Rlrg [%] 


SLRG [%] 


CC e+p 


1-8 1^ 




ARIADNE e+p 


0.918:? 


43 


MEPS e+p 


0.3 ±0.1 


39 


POMPYT e+p 


13.8 ±0.5 


53 


NC control sample (e+p) 


0.9 ±0.4 





Table 8: The fraction of large rapidity gap (LRG) events in the e+p CC sample, in the 
ARIADNE, MEPS and pompyt CC MC sets, and in the NC control sample. -Rlrg denotes 
the fraction of events with r] max < 2.5 and cos 6h < 0.75 (LRG events) in the total samples 
passing the CC selection cuts. For the NC control sample, -Rlrg is calculated from the 
sums of cross-section weights as described in section 3.3. The selection efficiency £lrg 
is defined as the ratio of the number of events which pass the CC cuts, fulfil the LRG 
requirement and have Q 2 COT > 200 GeV 2 to the number of events having r] max < 2.5 and 
Q 2 > 200 GeV 2 at the generator level. 
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Figure 1: Distribution of events in the (x, Q 2 ) plane. The filled (open) circles indicate 
the selected e + p (e~p) CC events in a sample with integrated luminosity of 2.93 pb -1 
(0.82 pb _1 ). The horizontal lines of the grid show the Q 2 , the vertival lines the x bins 
used for the analysis. The solid diagonal line indicates the kinematic limit y — 1. 
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Figure 2: Differential CC cross section versus Q 2 . In a), the filled (open) dots and the solid 
(dashed) curve represent the measured values and the Standard Model (SM) prediction of 
da/dQ 2 for e + p (e~p) collisions. Plot b) shows {da e + p ^ V x/dQ 2 )/{da e - p ^ v x/dQ 2 )- The 
error bars indicate the statistical and systematic uncertainties combined in quadrature. 
The horizontal position of the points is given by the generator level Q 2 average of the MC 
events in each bin. 
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Figure 3: Differential CC cross section versus x for Q 2 > 200 GeV 2 . In a), the 
filled (open) dots and the solid (dashed) curve represent the measured values and the 
Standard Model (SM) prediction of da/dx for e + p (e~p) collisions. Plot b) shows 
(da e + p ^ T ;x/dx)/(da e - p _ >u x/dx). The error bars indicate the statistical and systematic 
uncertainties combined in quadrature. The horizontal position of the points is given by 
the generator level x average of the MC events in each bin. 
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Figure 4: Differential CC cross section versus y for Q 2 > 200 GeV 2 . In a), the 
filled (open) dots and the solid (dashed) curve represent the measured values and the 
Standard Model (SM) prediction of da/dy for e + p (e~p) collisions. Plot b) shows 
{da e + p ^ v xl dy)/ {da e - p ^ vX / dy). The error bars indicate the statistical and systematic 
uncertainties combined in quadrature. The horizontal position of the points is given by 
the generator level y average of the MC events in each bin. 
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Figure 5: Differential CC cross section versus y for Q 2 > 200 GeV 2 , a) for x < 0.1 and b) 
for x > 0.1. The filled (open) dots and the solid (dashed) curves represent the measured 
values and the Standard Model (SM) prediction of da/dy for e + p (e~p) collisions. The 
error bars indicate the statistical and systematic uncertainties combined in quadrature. 
The horizontal position of the points is given by the generator level y average of the MC 
events in each bin. 
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Figure 6: Plots a) and b) show the distributions of transverse momenta and pseudora- 
pidities of all jets found in the e + p sample. Plots c) and d) give the transverse energy 
profiles in 77 and for jets in events with exactly one jet. Positive A77 values correspond 
to angles between the jet axis and the forward direction. The points denote the data and 
the solid (dashed) histograms represent the ARIADNE (meps) Monte Carlo predictions 
including detector effects. In plots a) and b) the MC predictions have been normalized 
to the number of jets observed in the data. 
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Figure 7: Plots a) and b) show the distributions of transverse momenta and pseudora- 
pidities of all jets found in the e~p sample. Plots c) and d) give the transverse energy 
profiles in rj and for jets in events with exactly one jet. Positive Ai] values correspond 
to angles between the jet axis and the forward direction. The points denote the data and 
the solid (dashed) histograms represent the ARIADNE (meps) Monte Carlo predictions 
including detector effects. In plots a) and b) the MC predictions have been normalized 
to the number of jets observed in the data. 
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Figure 8: Scatter plot of cos9 H vs. r] max containing all events in both the e + p and the e~p 
samples (see section 5.3 for the definitions of Oh and i] mSuX ). The region cos Oh < 0.75 and 
^max < 2.5, which corresponds to large rapidity gap events, contains a single event. 
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Figure 9: A display of the candidate large rapidity gap event. The left part shows the 
ZEUS inner tracking system and the calorimeter. The proton beam enters from the right; 
correspondingly, FCAL, BCAL and RCAL are displayed from left to right. The filled 
rectangles in the calorimeter denote energy deposits which are above the noise thresholds 
described in the text (cf. section 3.1) and in addition are above 100 MeV if they are 
isolated. The ZR-view shows a projection of tracks and energy deposits along constant 
values of Z and R = (X 2 + Y 2 ) 1 ^ 2 , where the two hemispheres are separated by a plane 
perpendicular to the event sphericity axis (20° < <fi < 200° corresponding to the upper 
half of the display). In the right part of the figure the ^-^-distribution of the calorimeter 
transverse energy is presented. The FCAL cluster at rj = 2.2 has 0.25 GeV of transverse 
momentum. The kinematic variables for this event are reconstructed to be f t —lA ±2 GeV, 
Q 2 =300 ± 70 GeV 2 , £=0.0093 ± 0.0015 and y=0.35 ± 0.10. 
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